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ABSTRACT 

We present high spatial resolution 35 parsec) 5 — 38^m spectra of the central region of M82, taken 
with the Spitzer Infrared Spectrograph. From these spectra we determined the fluxes and equivalent 
widths of key diagnostic features, such as the [NeII]12.8/im, [NeIII]15.5^m, and II2 S(l)17.03^m lines, 
and the broad mid-IR polycyclic aromatic hydrocarbon (PAH) emission features in six representative 
regions and analysed the spatial distribution of these lines and their ratios across the central region. We 
find a good correlation of the dust extinction with the CO 1-0 emission. The PAH emission follows 
closely the ionization structure along the galactic disk. The observed variations of the diagnostic 
PAH ratios across M82 can be explained by extinction effects, within systematic uncertainties. The 
16 — 18/im PAH complex is very prominent, and its equivalent width is enhanced outwards from the 
galactic plane. We interpret this as a consequence of the variation of the UV radiation field. The 
EWs of the 11.3/im PAH feature and the H2 S(l) line correlate closely, and we conclude that shocks 
in the outflow regions have no measurable influence on the H2 emission. The [Nelll] / [Nell] ratio is 
on average low at ~0.18, and shows little variations across the plane, indicating that the dominant 
stellar population is evolved (5-6 Myr) and well distributed. There is a slight increase of the ratio 
with distance from the galactic plane of M82 which we attribute to a decrease in gas density. Our 
observations indicate that the star formation rate has decreased significantly in the last 5 Myr. The 
quantities of dust and molecular gas in the central area of the galaxy argue against starvation and for 
negative feedback processes, observable through the strong extra-planar outflows. 

Subject headings: galaxies: starburst — galaxies: individual (M82) — infrared: galaxies 



\r' 1. INTRODUCTION 

> ■ 

(3 ' M82 ( NGC 3034) is an irregular galaxy located at 3.3 
^ \ Mpc (Freedman fc M adorc 1988) in the M81 group. It 
' is the closest starburst galaxy, seen nearly edge-on, with 
■ an inclination angle of about 80deg. At infrared wave- 
' , lengths it is the brightest galaxy on the sky, with a tota l 
^ ' infrared luminosity of 3.8 x W°Lq (jColbert et al.lll999D . 
J — Most of its luminosity originates from the inner 500 pc 
, hosting intense starburst activity pres umably triggere d 
by a tidal interaction with M81 (e.g. lYun et ahl I1993D . 
Evidence for this interaction comes from the HI streams 
K> ■ which connect M81 to all three (M82, NGC 3077, and 
^ ' NGC 2976) peculiar members of the inner M81 group 
C3 ' ()Appleton et al.ll 19811 : IVun et al.lll994f ). Recent deep op- 
tical images also revealed star s associated wit h the HI 
bridge between M81 and M82 (|Sun et al.ll2'005h . 

At the distance of M82, 1" corresponds to 15 parsec, 
which allows spatially resolved studies of the starburst 
region. Evidence of a stellar bar ~ 1 kpc long is 
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shown by near-inf rared studies fe.g. iTelesco et aLlll991l : 
iLarkin et al.lll99^ . mid-infr ared [NeII]12.8 um and mil- 
lii netric CO emission s tudie s (iLo et al.lll987l). Ac c ordin g 
to ILarkin et all (|1994( l and lAchtermann fc Lacvl (|1995( l. 
there is a rotating ring of ionized gas at a radius of 
~ 85 pc, and on the inner side of a ring of molecular 
gas at ~ 21 pc. Two possi ble s piral arms were als o 
identified bv lShen fc Lol (Il99l and lMavva eraD ((2005h . 
at radii of ~ 125 pc and ~ 400 pc. The starburst 
of M82 drives a bipolar mass outflow out to several 
kiloparsecs perpendicular to the plane of the galaxy, 
espec i ally evident in X-ray and Ha (iBregman et al 



1^ iShopbell fc Bland-Hawthoral 119981: iLe^ert et al 
1999t iCappj et al.l 119991 : [Strickland etTaH [200I . 
Dust has also been detected in the outflow region 
(lAlton et al.l 119991: iHeckman et al.l 120001 : iHoopes et all 
120051 lEngelbracht et all l2006l ). The star forming 
regions of M82 are predominantly clustered in the 
volume enclosed by the molecular gas ring, indicated 
by the HII region tracers, like the [NeII]12.8/im line 
and the mid - and far-infrared c ontin u um emission 
(iTelesco et al.l 119911: IWalter et al l 120021: ILarkin et al. 
19941 lAchtermann fc Lacvl 119951 : iLipscv fc Plavchan 



20041 ). Near-infrared hydrogen recombination lines also 
arise in these regions, but they are a more ambiguous 
tracer of star formation, as they can also be excited by 
shocks, though these are unlikely to dominate. 

Near-infrared integral field spe ctroscopy and ISO-SWS 
mid-in frared spectroscopy by iForster Schreiber et al.l 
(|200lD allowed a detaile d modelling of starburst actiy - 
ity in the central region (jForster Schreiber et al.l l2003b'). 
These models are consistent with the occurrence of star- 
burst activity in two successive episodes, about 10 and 
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Fig. 1. — Left: Overlay of SL (red) and SH (blue) coverages 
and selected low-res extraction regions on an IRAC image. 
Regions A and B, represented in yellow, are regions where LL (14 
- 35 ^im) spectra were extracted. Right: Zoom-in of the IRAC 
8fim image with an overlay of the SH map area in blue and the 
selected regions (in green) from where the SH-I-SL spectra were 
ex tracted. The ISO-SWS apertu re used in the 12 — 27/jm range 
bv lForster Schreiber et al.l I I2001I ) is also overlayed in yellow on the 
image. Both figures are in logarithmic scaling. 

5 Myr ago, each lasting a few milli on years . How ever, 
the spatial studies by Forster Schrei ber et al.l (l2001[ ) and 
others covered only near-infrared wavelengths. The 
large aperture of the ISO-SWS provided a continuous 
2.4— 45/im spectrum but covered the whole central region 
of M82. ISOCAM-CVF data ijForster Schreiber et all 
[2003a) provided better spatial resolution but the spec- 
tra only had a spatial resolution of i? 40 and were 
s hortward of 15/im. 

lEngelbracht et all (|2006l ) published Spitzcr-IRS low- 
resolution spectra of a 1 arcmin wide strip along the mi- 
nor axis of M82, intersecting the disk at the eastern side. 
The spectra, taken as part of the SINGS Legacy project, 
were combined with Spitzer 8 and 24 fim images, and 
show that the emission by polycyclic aromatic hydro- 
carbons (PAH) and H2 molecules extends f ar out from 
the di sk (to 6 kpc) in both directions. Engel bracht et all 
(|2006) suggest that the filamentary aromatic-dominated 
emission represents dust either expelled from the galaxy 
as a result of a powerful nuclear superwind, or that dust 
is in the halo being lit up by the starburst, perhaps co- 
existing with the extensive warm H2 molecules. They 
suggest that this halo dust is probably a leftover from 
the interaction with M81. 

In this paper we present mid-IR spectral maps at un- 
surpassed sensitivity and spatial resolution of the central 
~ 0.5 kpc^ of M82, covering the main contributors to 
the bolometric luminosity of the galaxy. Our goal is to 
provide a spatial and spectrally detailed description of 
the physical conditions within the central ~ 500 pc of 
M82, to help us to give an insight on the evolution of the 
starburst activity in this region. This involves the study 
of the distribution of the radiation field, gas density, and 
the physical properties of PAHs. Of particular interest 
are the spatial variations of the fine-structure lines, the 
excitation of the molecular hydrogen, and the distribu- 
tion of the PAH molecules. In section 2 we describe the 
observations and data reduction, in section 3 we present 
the data, and in section 4 we discuss the scientific results, 
followed by our conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were made with the Infrared S pec- 
trograph (IRS)^ spectrometer (jHouck et al.l I2004D on 
board the Spitzer Space Telescope, under the IRS guar- 

* The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation funded by NASA through the 
Jet Propulsion Laborator and the Ames Research Center 



anteed time observing program. The data were taken 
on June 6th, 2005 using the IRS "mapping mode" in all 
four modules: Short-High (SH; 10 — 19/J,m), Long-High 
(LH; 14 - 38/im), provide R ~ 600, while Short-Low 
(SL; 5 — 14/im) and Long-Low (LL; 14 — 38/im) give 
i? ^ 60 — 130. Each of the SL and LL modules are 
further divided into two subslits, which correspond to 
diffraction orders: SLl (7.5 - 14/im), SL2 (5 - 7.5//m), 
LLl (20 - 38/im), and LL2 (14 - 20/im). The SH map 
consists of 30 pointings with 4 cycles each, and each sub- 
sequent pointing is offset by half a slit width parallel to 
the slit and about one third of the slit length along the 
slit. The SH map covers an area of 28" x 23". The LH 
data consists of 12 pointings, with 5 cycles each, covering 
an area of 38" x 33". The offsets are equivalent to the 
SH map. The SL data consists of 120 pointings, with 2 
cycles each, covering an area over the M82 central region 
of 55". 5 X 57". The LL data consists of 22 pointings with 
2 cycles each. Both SL and LL maps follow the same off- 
seting scheme as SH. Because of the high brightness of 
M82 it was unavoidable that the LLl data became satu- 
rated near the center of M82, but the LL2 data are still 
usable. Fig. [T] (left) shows the areas covered by the IRS 
SL, LH, and SH maps overlayed on the IRAC 8/im im- 
age from [EngelbrachtIeF|n] ([200^. The total integration 
times range from 12 (for SL exposures) to 31 sec (for LH 
exposures). Due to issues concerning the extraction, the 
LH spectra were not used in this analysis. The boxes 
'A' and 'B' are the regions where the complete low-res 
(SL-I-LL) spectra were extracted. 

The data were processed with version 13.2 of the 
Spitzer reduction pipeline (version 14 for LL). Observa- 
tions taken at each po sition were combin ed into spectral 
cubes using CUBISM (jSmith et al.ll2007t ). an IDL-based 
software package designed to combine spectral mapping 
datasets into 3D spectral cubes. Bad pixels in the ba- 
sic calibrated data (BCD) spectra were manually flagged 
and then automatically discarded when rebuilding the 
cube. Spectra from off-source positions 1 kpc to the NE 
of the nucleus were used to subtract the background from 
the low-resolution spectra. For the high-resolution spec- 
tra we did not subtract a background since there was no 
suitable "sky" spectrum available and the high source 
fluxes strongly dominate any background emission. The 
spect ral analysis wa s done using SMART (iHigdon et all 
[200I and PAHFIT (|Smith et al.ll2006a|) . 

3. ANALYSIS 

On the basis of the SH spectral map coverage we de- 
fined six sub-regions for which we extracted the spectra 
from SL and SH with CUBISM. The location of these 
regions is shown in Fig . [1] (right) ove rl ayed on the IRAC 
8/im image from lEngelbracht et all (|2006| ). The ISO 
aperture is also shown to illustrate the increase in spatial 
resolution obtained with Spitzer-TRS. Regions 2 and 3 
coincide with the peaks of the [Nell] emission, regions 1 
and 4 are offset to both sides of the galactic disk, the 
slightly larger region C covers the nominal center of the 
galaxy, and the last region corresponds to the entire area 
mapped with the SH spectrograph. The exact coordi- 
nates and sizes are listed in Tabled 

Fig. [2] shows the SL (5 - 14/im), and SH (10 - 19/im) 
spectra extracted for these six regions. The noise 
is negligible and the spectra exhibit the classical fea- 
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TABLE 1 

Characteristics of the Selected Regions 





Region 1 


Region 2 


Region 3 


Region 4 


Center 


Total 


Region A 


Region B 


RA (J2000) 
Dec {J2000) 
Size 


9:55:52.42 
+69:40:32.1 
6".8 X 6".8 


9:55:50.63 
+69:40:45.6 
6". 8 X 6". 8 


9:55:52.74 
+69:40:48.5 
6".8 X 6".8 


9:55:50.22 
+69:40:58.60 
6". 8 X 6". 8 


9:55:51.19 
+69:40:46.8 
11".3 X 11".3 


9:55:51.32 
+69:40:45.30 
24".9 X 24".9 


9:55:52.17 
+69:41:12.5 
14".8 X 20".3 


9:55:52.68 
+69:40:22.2 
14".8 X 20".3 



tures of starburst galaxies (|Brandl et alJ l2006l ). such 
as strong emission features of the PAHs, fine-structure 
hues, and emission from molecular hydrogen, in addi- 
tion to the broad silicate absorption features around 
9.7/im and 18^m, with an underlying continuum of 
very small grain emission (VSGs). The 5-38 /zm wave- 
length range contains many important diagnostic lines, 
such as [ArII]6.99/im, [ArIII]8.99/im, [SIV]10.51/im, 
[NeII]12.81^m, [NeIII]15.56Aim, [SIII]18.7lAim. AUofthe 
detected features are labelled in Fig. [2] The fluxes and 
ratios of the most relevant fine-structure and H2 emission 
lines are listed in Table O These lines were measured us- 
ing Gaussian fits to the line and linear fits to the local 
continuum. The [Nell] line was measured after subtract- 
ing the f2.6/im PAH feature. [SIV] at 10.5/im is hard 
to detect on regions 1 and 4 due to low S/N ratio, and 
[Arlll], at 8.99/im is too close to the 8.6/im PAH feature 
to be detected in a low resolution spectrum. Due to sat- 
uration of the LLl module and extraction issues of the 
LH spectra, we could not measure the [SHI] 33.6/im line. 

Numerous PAH emission features are easily detectable 
in our spectra. The fluxes and equivalent widths 
(EWs) of the strongest features at 6.2/im, 7.7/im, 8.6/im, 
11.3/im, 12.6/im, 14.2/im, and the 16 — 18/im complex 
are listed in Table [3l Their values were derived using 
PAHFIT (jSmith et al.ll2006a[ ). an IDL tool that decom- 
poses low-resolution spectra of PAH emission sources us- 
ing a physically motivated model. This model includes 
starlight, thermal dust continuum, resolved dust features 
and feature blends, prominent emission lines, and dust 
extinction. In our case, we merged SL and SH spectra 
in order to have one spectrum for each region with the 
widest wavelength coverage possible. Weaker PAH fea- 
tures at 5.2 /im, 5.6 /zm, 12.0 /im, and 13.55 fim are 
present in the spectra in Fig. [21 but are not further 
utilised in this paper. In Fig. [3] we present an exam- 
ple of a fit to a combined SL-t-SH spectrum of region 2. 
In the overall spectrum (green) we distinguish the con- 
tinuum component (red), PAH component (blue), ionic 
lines (purple), and dust absorption (dashed line). 

Silicate absorption affects mainly the wavelength inter- 
val from 7 — 12/im. In PAHFIT we have the option to 
include or exclude silicate absorption, and the difference 
between these two cases can be up to 50% in flux of the 
11.3/im feature. In section 4.3 we discuss the methods to 
estimate the magnitude and distribution of extinction. 
Also there is a significant difference in the 11.3/im flux 
between the SH and SL spectra, which can also reach 
50%. This difference can be attributed to a poor fit to 
the silicate absorption feature in SH, as this module only 
covers wavelengths > 10/im. In this paper we use the 
SH+SL measurements for this reason and to include the 
17/im complex. 

We also extracted combined SL-I-LL spectra (5 — 38/im) 
from two 15" x 20" regions, A and B (see Fig. 1 left), 
at a distance of 200 pc above and below the galactic 



plane of M82, where saturation did not compromise the 
LL measurements. Their positions are listed in Table 1. 
The spectra from regions A and B were virtually identi- 
cal in shape, and an average of these spectra is shown in 
Fig.m The spectrum is clearly dominated by strong PAH 
emission features and a steeply rising conti nuum, charac- 
teristic of classical starburst galaxies (e.g. iBrandl et al.l 
_2006). The wiggles observed longwards 20/im are fringes 
that originate from interferences in the detector substrate 
material. At the low spatial resolution presented by the 
LL modules, we cannot see any unusual features in the 
SL+LL spectra of regions A and B. Also, as the LLl slit 
is saturated at the central region of M82, we cannot ex- 
tract any full low-res spectra of this region, connecting it 
to regions A and B. For these reasons, we will not make 
further analysis of LLl module spectra in this paper. 

4. RESULTS AND DISCUSSION 

The main gain of our observations over previous work 
on the starburst in M82 is the availability of spatially 
resolved mid-IR spectroscopy of the central region. De- 
spite their overall similarity the spectra show distinct 
variations in the relative strengths of the spectral fea- 
tures, in particular in the neon fine structure lines and 
PAH features. These variations and their physical causes 
will be described in the following sections. 

4.1. The Morphology of the Starburst Region 

The discovery of a series of compact radio supernova 
remnants al ong the galactic plane of M82, ext ending 
over 600 pc (jKronberg et al.lll985HMuxlow et al] [r994). 
is an indication of very recent and presumably ongoing 
star formation. The detection of the ionic high excita- 
tion lines [Nelll] and [SIV] confirms the presence of very 
young massive stars in M82. Ratios using ionic lines of 
the same species and different ionization potentials such 
as [NeIII]/[NeII], [SIV]/[SIII], and [ArIII]/[ArII] are a 
useful measure of the hardness and intensity of the radi- 
ation field and radiation density, and are therefore sen- 
sitive to the presence of young massive stars. 

Our measurements of the [NeIII]/[NeII] ratio are shown 
in Table [2l for each of the six selected regions. We 
find 0.13 < [NeIII]/[NcII] < 0.21, with a median of 
[NeIII]/[NeII] — 0.18, which is consistent with the spa- 
tially integra ted ratio of [Nelll] / [Nel l ] = 16 ± 0.04 de- 
termined bv lF5rster Schreiber et al] ( 20011 ). These val- 
ues are 30% lower than the average valu e of 0.26 for the 
ISO- SWS sample of starburst galaxies (jThornlev et al.l 
[2000l) . and more than an order of magnitude below 
the value of 8.5 found in the center of NGC 5253 
(jBeirao et al.ll2006f ). a low metallicity starburst galaxy 
at about the same distance as M82. 

Other ratios such as [SIV]/[SIII] and [ArIII]/[ArII] 
could be used to confirm the results on the [NeIII]/[NcII]. 
However, as mentioned in Chapter 3, our measurements 
of [SIV] in regions 1 and 4 have large errors from noise. 
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Fig. 2. — SL (left) and SH (right) spectra of the six selected regions within the central starburst region of M82. The "Total" region 
corresponds to the total SH map in Fig. [T] 
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Wavelength (|j.m) 

Fig. 3. — Decomposition of SL+SH spectrum of Region 2. Red 
solid lines represent the thermal dust continuum components, the 
thick gray line the total continuum, blue lines are dust features, 
while the violet peaks are atomic and molecular spectral lines. 
The dotted black line indicates the fully mixed extinction which 
affects all components, with axis at right. The solid green line is 
the full fitted model, plotted on the observed flux intensities and 
uncertainties. 

and [Arlll], at 8.99/im is too close to the 8.6fim PAH fea- 
ture to be measured accurately from the low resolution 
spectrum, precluding the use of these ratios. 

The Spitzer/IRS spectral maps allow a study of the 
spatial variation of the radiation field, based on the 
[NeIII]/[NeII] ratio. Due to the presence of several lumi- 
nous star clusters in the central region, one might expect 
strong variations of [NeIII]/[NeII] between regions. In 
Fig. [5] we present spectral maps of the two strongest neon 
emission lines, and a ratio map, overlayed by th e 12/xm 
continuum contours from Acht ermann fc Lacvl (|1995[ ). 
For both line maps, CUBISM was used to subtract 
a fitted continuum map from a total line-|-continuum 
map at the same wavelength on a pixel-by-pixel ba- 
sis. The [Ne II] emission in the upper map shows two 
peaks to either side of the nucleus, a weaker and a 
st ronger peak, which correspo nd to the E and W peaks 
in I Achtermann fc Lacvl (|1995f ) respectively. After a close 
inspection, we identify the strong W peak in our maps 
with the two continuum emission peaks. The E and 
W [Nell] emission peaks are identified as a "ring" of 
ionized gas in Achtermann fc Lacv (199_^). The mor- 
pholog y of the [Nell ] emis sion foll ows the Br7 emis- 
sion inlSatvapal et al.l (|1997l ) and iForster Schreiber et all 
(|200ll) . The [Nelll] map also reveals two emission peaks 
at the same positions, but with the eastern peak brighter 
relative to the western peak. 

The lower map in Fig. [5] presents the [NeIII]/[NeII] ra- 
tio. The ratio varies from 0.08 — 0.27 throughout the 
map. The lower value corresponds to the location of the 
westernmost clusters and there is a significant increase in 
the ratio further out from the galactic plane of M82, from 
0.15 to 0.27. Statistical la errors are ~ 15% and arise 
mostly from baseline determination errors. The grey line 
repre sents the direction of the X- ray outflow observed in 
M82 (jStrickland et al.lll997l . l2004D . It originates from the 
nucleus (marked with a plus sign) and is perpendicular 
to the plane of the galaxy. The peak of X-ray emission 
is offset from the outflow axis by 30 pc. Although the 
gradient in the [NeIII]/[NeII] ratio does follow the out- 
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Fig. 4. — Average of the 5 — 38/xm low-resolution IRS spectra 
of the regions A and B, located approximately 200 pc above and 
below the galactic plane of M82. 

flow axis, it appears offset by 5 arcsec to the east, and 
associated with the eastern cluster. 

The overall morphology of the ionizing radiation in 
the region as revealed by the spectral maps seems to 
be in agreement with ground observation s of the contin- 
uum and [Nell] emission by Achtermann fc Lacvl |i995). 
However, the [NeIII]/[NeII] ratio varies only by a fac- 
tor of three throughout the region and these variations 
do not correspond to the position of the emission peaks. 
An increase of the [NeIII]/[NeII] ratio is observed fur- 
ther out of the galactic plane, but due to the limited 
spatial coverage of the map, we cannot determine if this 
increase is connected to the outflows. The origins of the 
observed emission morphology and the variation of the 
[NeIII]/[NeII] ratio are discussed in the following subsec- 
tion. 

4.2. Origins of the Variation of the Radiation Field 

The [NeIII]/[NeII] ratio measures the hardness of the 
radiation field, which is a function of stellar age and 
metallicity, and often parameterized by the effective tem- 
perature Te//, and the radiation intensity as measured 
by the ionization parameter U^. Given the strength of 
the [Nell] and [Nelll] lines in M82, and their proximity 
in wavelength which minimises the effects of extinction, 
[Nelll] / [Nell] is the most reliable measure of the hard- 
ness of the radiation field. With the help of starburst 
models, it is possible to use the [Nelll] / [Nell] ratio to 
estimate the ages of the massive clusters in the region. 

Observations with IS O/SWS have been used previ- 
ously for this purpose. iForster Schreiber et al.l ()2001l ) 
determined a spatially integrated ratio of [Nelll] / [Nell] 
= 0.16 for the inner 500 pc of M82. Using the 
photoionization code CLOUD Y and solar m etallicity 
stellar atmosphere models b y iPauldrach et al.l (|1998l ). 
iForster Schreiber et al.l (|2001l ) modelled the variations of 
line ratios with Teff. Adopting an electron density — 
300 cm~^ and an ionization p arameter \ogU = —2.3, 
IForster Schreiber etiill (|2003bh found for the ISO value 

^ U is defined as U = ^^p^^ — where Q is the production rate 
of ionizing photons from the stars, R is the distance between the 
ionizing cluster and the illuminated gas cloud, njj is the hydrogen 
number density of the gas, and c is the speed of light. 
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Fig. 5. — Spectral maps in the [Nell] (top), [Nelll] (cen- 
ter), and [NeIII]/[NeII] (bottom) lines from the IRS SH mod- 
ule, with c ontour ov e rlays of the 12/^m continuum emission from 
lAchtcrman n &: Lacvl 119951) . The regions shown arc ~ 26" X 26". 
At the bottom map, the line is the direction of the of the X-ray 
outflow, which is perpendicular to the plane of the galaxy. The 
cross represents the nucleus. 

for [NeIII]/[NeII] an effective temperature of 37400±400 
K and a burst age of 4 -6 Myr. Other observed ratios were 
also modelled, such as [ArIII]/[ArII] and [SIV]/[SIII], 
giving similar T^^ff (within uncertainties). Independent 
estimati ons of cluster ages were done at longer wave- 
lengths. IColbert et al.l fl999i ) analysed far-infrared spec- 
tra from ISO/LWS and fitted line ratios to a combined 
HII region and PDR model. Their best fit model is an 
instantaneous star burst of 3 - 5 Myr old massive stars, 
in agreement with lF5rster Schreiber et al.l (|2003br ) . 

If the infrared emission peaks correspond to mas- 
sive clusters of stars, we can determine their ages from 
the measured [Nelll] / [Nell] ratio and compare them 
with t he above result s , usin g the photoio nization mod- 
els by iSniiders et al.l ()2007f ). As input, 'Sniid ers et al.l 
(|2007l) used massive cluster spectra modeled with Star- 



burst99, assuming a Salpeter IMF, M„ 



lOOM, 



0: 



Mlow = O.2M0, and a gas density of 100 cm"^. This 
value of gas density is a fact or of three lower than the 
iForster Schreiber et al.l (|2001h value, but this has a small 
effect on the [NeIII]/[NeII] ratio. Fig. [6] show the results 
for a range of ionization parameters. The selected re- 
gions are represented by horizontal lines. The value o f 
logU = -2.3 derived by iForster Schreiber et al.l (|200lD 
corresponds to a ionization front speed oiq = 1.6 x lO^cm 
s^^. For a typical value of q, the ages of the clusters in 
each region range from 5-6 Myr, in agreement with the 
previously determined burst ages. These clusters dom- 
inate the central region of M82 and may be similar to 
those obser ved further out in M 82 in the optical, studied 
in detail bv lSmith et al] (|2006b[ ) using HST/ACS, which 
were found to have an average age of 6.4 ± 0.5 Myr. 

As shown in Fig. [51 the [NeIII]/[NeII] ratio increases 
from 0.15 to 0.27 with increasing distance from the galac- 
tic plane of M82. This is counterintuitive, as one might 
expect a harder radiation field at the location of the most 
luminous regions along the plane. Away from the plane, 
a decrease of gas density, relative to the number of ioniz- 
ing photons, leads to an increase of the ionization param- 
eter, which then cause s an increase of t he [Ne III]/[NeII] 
ratio, as discussed in iThornlev et al.l ()2000[ ). Indeed, 
Fig. [6] shows that the [NeIII]/[NeII] scales with the ion- 
ization parameter. The variation we observe in the 
[Nelll] / [Nell] ratio implies a variation of a factor of five 
in the ionization parameter. This is equivalent to say- 
ing that the gas density decreases five times faster than 
the radiation field which decreases as ~ with R be- 
ing the distance to the ionization source. Shocks could 
also contribute to the increase of the [NeIII]/[NeII] ra- 
tio in the outflow region, but in Section 4.5 we show 
that to be minimal. This hypothesis could be tested us- 
ing the [SIII]18.6/im/[SIII]33.6/im. Unfortunately, due 
to the problems reported in Sec. 2, we could not derive 
an accurate flux for the [SIII]33.6/im line in both LL and 
LH spectra. 

It is important to emphasize that, even at higher an- 
gular resolution, the [NeIII]/[NeII] ratio in M82 remains 
quite low for an active starburst. We would have ex- 
pected a larger variation with higher ratios locally cor- 
responding to younger clusters and lower elsewhere. A 
con iparison with th e ISO- SWS sample of starburst galax- 
ies (jThornlev et al.l (120001 )) shows that it is actually lower 
than most starbursts, despite being closer and better 
resolved. The low [NeIII]/[NeII] ratio could be caused 
by an aged stellar population where the starburst activ- 
ity ceased more than half a dozen Myrs ago - although 
this possibility seems unlikely given the large amounts 
of molecular gas still present at the center of M82. An 
edge-on view of the galaxy could also contribute to the 
low variation of [NeIII]/[NeII] ratio. 

As a comparison, we examined the 4". 5 x 4". 5 area 
with the highest [Nelll] flux (Fig. O, comparing the 
measured [Nell] and [Nelll] luminosities with the mod- 
els to determine the enclosed stellar mass. A single su- 
per star cluster of 5 Myr would have a cluster mass of 
1O^M0 which is twice the mass o f the super star cluster 
in NGC 5253 ([Turner et al.ll2003f) . A single cluster in a 
4". 5 X 4". 5 area would correspond to a cluster number 
density of 200/kpc^, which is comparable to the clus- 
ter density found in the fossil starburst region of M82 by 
Ide Griis et all (pO?)l . 
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While the average age of the starburst population in 
M82 appears to be ~5 Myr, ongoing star formation (<1 
Myr) could possibly be obscured by recent contributions 
of older stellar populations (> 5 Myr) to the neon ratio. 
Considering this possibility, we explore the region with 
the highest [Nelll] flux, a strong indicator of the pres- 
ence of O stars, as an illustrative case to set an upper 
limit on the ongoing star formation in M82. This region 
has a [NeIII]/[NcII] of ~ 0.10, with a total [Nelll] lumi- 
nosity of 1.41 X 1O°L0. Assuming ^ 1.6 x lO^cm s~^, a 
typical young cluster of 1 Myr has [NeIII]/[NeII]~ 5 and 
can contribute ~50% to the total [Nelll] luminosity or 
~2.5% per mass relative to the older 5—6 Myr old pop- 
ulation. This means that the 1 Myr old population only 
emits ~5% of the total [Nell] -t- [Nelll] luminosity of the 
older 5—6 Myr old population. Even as an upper limit, 
this value indicates that the activity of the starburst has 
substantially declined relative to the high star formation 
rate that existed 5 Myr ago. 

Despite the reduced starburst activity, the presence 
of CO emission all over the central region (Fig. [7]) 
shows that the starburst in M82 still has a larg e gas 
reservoir, as pointed out bv iThornlev et all ()2000l ) and 
iForster Schreiber et all (|2003bf ). It appears therefore un- 
likely that the starburst activity ceased because of lack 
of fuel. On the other hand, negative feedback mainly 
through strong stellar winds and supernovae explosions 
can play an important r ole in determining the sta r forma - 
tion rate in starbursts. iForster Schreiber et al.l ()2003bl ) 
calculated a feedback timescale of 1 — 10 Myr, which is 
in good agreement with the age of the older super star 
clusters in the central region of M82. 

In summary we find that [NeIII]/[NeII] ratio is low 
on average, and increases with distance from the galac- 
tic plane of M82. The increase can be explained by an 
increase of the ionization parameter through a drop in 
gas density. The low [Nelll] / [Nell] ratio indicates that 
the dominant population consists of older clusters (> 5 
Myrs). We cannot rule out the presence of younger clus- 
ters in the central region, but at a much reduced rate 
(< 5%) of star formation compared to previous epochs. 
The large amount of molecular gas still present in the 
central region argues against a starvation of starburst 
activity, due to lack of gas. Instead, we believe that neg- 
ative feedback processes are responsible for the observed 
decline in the star formation rate. 

4.3. Extinction 

One of the most noticeable features in the spectra in 
Fig. [2] is the absorption feature at 9 — 11/xm, caused by 
silicate grains. This f eature can be extre mely deep espe- 
cially in ULIRGS (e.g lSpoon et al.ll2007[ ). and affects the 
observed fluxes of spectral lines and features in this re- 
gion. The strength of this feature is characterized by the 
optical depth at 9.8/^m (rg g) assuming a simple geomet- 
rical dust distribution. For our analysis it is important 
to estimate the intensity of this feature and its spatial 
variation in order to investigate its influence on the PAH 
strengths, which we discuss in the following subsections. 

To study the distribution of rg.g in the SL region, we 
selected an area of 20 x 12 pixel in the center of the SLl 
map. Within that area we extracted spectra from 60 
spatial resolution elements (2x2 pixel each). 

The apparent optical depth rg.g is then estimated from 
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Fig. 6. — Effect of the cluster age on the [NeIII]/[NeII] ratio. 
The data for our selected regions in M82 is shown by the hor- 
izontal lines. The model curves are computed for a cluster of 
10'' Mq, assuming Salpeter IMF, Mup = lOOAf0, and = 
0.2A/q. Each curve represents a different ionization parameter, 
and the solid curve is the one that approaches the value found by 
[Forster Schreiber et al. (200]]), logU = -2.3. The horizontal lines 
indicate the [NeIII]/[NeII] for each selected region, from region 1 to 
region 4. The first dip in [NeIII]/[NeII] represents the ageing of the 
stellar population, after which WR stars are produced, increasing 
the [NeIII]/[NeII] ratio. The dip at 6 Myr occurs as most massive 
stars die through supernova explosions. 

the ratio of the interpolated continuum to the observed 
flux at 9.8fim. For the central region of M82, rg.g ranges 
from 0.3 - 3.1, assuming a foreground screen attenuation. 
We built an extinction map of 10 x 6 resolution elements 
from the simple fit method, which is shown in Fig. [71 The 
qualitative distribution of rg.g observed in this figure is 
very similar to the qualitative results from the PAHFIT 
fitting method. There is a good correlation between the 
rg.g distribution and the CO 1-0 emission, indicating that 
dust and molecular gas coincide in this region. The en- 
hanced rg g in the northwest region of the map indicates 
an increase of silicate dust above the galactic plane. 

For an independent, and possibly more accurate, es- 
timate of rg.g, we used PAHFIT on 15 SL+LLl spec- 
tra (5 — 20/zm), each corresponding to 4 resolution ele- 
ments fitted only with SL. With the PAHFIT method, 
rg g ranges from - 2.52, with a median is 1.34. 

The combined SL-I-LL spectrum is necessary to better 
constrain the parameters in PAHFIT. Fitting only SL 
spectra with PAHFIT can result in large errors in the 
calculation of rg.g. Unfortunately, due to saturation at 
wavelengths longward of 20/im, we have an insufficient 
spatial coverage of LL data in the central region. In 
addition, the LL slits are ^ 10". wide, providing low 
spatial resolution. 

For the higher sp atial resolution, we used the simple 
method described bv lSpoon et al.l (|2007D to estimate rg.g. 
In this method we approximate the mid-IR continuum at 
9.8/im by a power law fit to the flux pivots at 5.5/im and 
14.5//m, avoiding the PAH emission features. 

The values derived from both SL-f LL and simple fit 
methods agree qualitatively well, but in regions where 
rg.g < 1, the difference between the two methods is 
greater than a factor of two. Using a mixed attenua - 
tion law with ^^/rg.g = 16.6 (Rieke &: Lebofskylll985( l. 
the range of rg.g corresponds to < A„ < 41.8 for the 
SL+LL method. The values of rg.g from the power-law 
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Fig. 7. — Map of rg g, determined from a simple continuum fit 
to the SL s pectra, with overlay of CO 1-0 emission contours from 
IShen &: Ld ||T999!). The crosses indicate the [Nell] emission peaks. 
The image was rotated by 55 degrees and interpolated. 

interpolation method gi ve 5.0 < Ay < 51.5. Th e se va l- 
ues are in agreement with lForster Schreiber et al.l (|200in , 
which derive 23 < ^„ < 45 for their selected regions in 
M82, which cover an area closer to the infrared peaks. 

We conclude that both methods, besides their signif- 
icant uncertainties in the magnitude of rg.s, have con- 
sistently revealed significant variations in the amount of 
dust extinction across the central region. These varia- 
tions are strong enough to affect the following interpre- 
tation of PAH features. 

4.4. Variations of PAH Emission Features 

Polycyclic Aromatic Hydrocarbons (PAHs) are 
thought to be responsible for a series of broad emission 
features th at dominate the mid- infrared spectra of star- 
bursts fe.g. iPeeters et alll2Q04bD . They are observed in 
a diverse range of sources with their strongest emission 
originating in photodissociation regions (PDRs), the 
interfaces between HII regions and molecular clouds. 

The relative strength of the different PAH bands is 
expected to vary with the size and the ionization state 
of the PAH molecule (Draine & Li 2001, 2 00j). Obser- 
vations of Galact i c sources (e.g. IVerstraete et al.l 119961 : 
iJoblin et al]|1996l : IVermeii et al.ll2002D have shown that 
the relative strengths of individual PAH features depend 
upon the degree of ionization: C-C stretching modes at 
6.2/im and 7.7/im are predominantly emitted by PAH 
cations, while the C-H out-of-plane bendin g mode at 
lL3um arises mainly from neutral PAHs (jPraine fc Lil 
120011 ). Thus the ratios 6.2/11.3Aim and 7.7/11.3/zm 
may be use d as indicators of PAH ionization state. 
IJoblin et all ([1996) have found that the 8.6/11.3 fim ra- 
tio can also be linked to variations in the charge state of 
the emittin g PAHs, which is supported by la boratory ex- 
perime nts ()Hudgins fc Allamandolal[l995l ). ISmith et al.l 
()2006a ^ found band strength variations of factors of 2-5 
among normal galaxies. 

Studies of PAHs in M82 have been done pre- 
yiously using ISO. Observ ations with ISOCAM 
(|Forster Schreiber et al.l l2003a[ ) revealed a decrease in 
the 6.2/7.7/im ratio and an increase in the 8.6/11.3/im 
PAH ratio from the nucleus outwards along the galactic 
plane. The observations also showed a good spatial 
correlation of the 8.6/11.3/im ratio with the CO (1-0) 
emission. These ratio variations are attributed to real 
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Fig. 8.— IRS spectral map of the PAH ratio 6.2/11.3^(m, with 
overlay of CO 1-0 emission contours from Shen & Lo (1995). The 
map is in logarithmic scaling. This map was done using CUBISM 
maps of baseline-subtracted flux at 6.0 — 6.5/xm and 11 ~ 11.7/tm. 

differences in the variation of physical characteristics of 
PAHs across M82, specifically a higher degree of PAH 
ionization within the most intense starburst sites. 

We use CUBISM to build maps of ratios of PAH fea- 
tures that exist in the 6-14 /im SL spectra. The 7.7/xm 
feature is split between the two SL orders, and the 8.6/im 
feature is largely influenced by the 7.7/xm feature, making 
the local continuum determination difficult. For these 
reasons the 6.2/im feature map was chosen to be com- 
pared to the 11.3/im feature map. In Fig. [8] we present 
the 6.2/11.3/im PAH r atio map , overlayed with the CO 
(1-0) contours from Sh en fc Lol (j 19951 ). 

Obviously, the distribution of the 6.2/11.3/im PAH ra- 
tio correlates well with the molecular ring defined by the 
CO (1-0) map. The molecular emission forms two lobes 
~ 200 pc from the nucleus whereas the ionized emission 
is concentrated around the nucleus, in the regions where 
the [Nell] and [Nelll] emissions peak. 

However, the 6.2/11.3/im PAH ratio may be affected 
by extinction. Comparing the CO emission with the ex- 
tinction map in Fig. [7] and with the PAH ratio map we 
see that they correlate well. Hence, we investigate the 
effects of extinction on the observed PAH ratios in the 
next subsection. 

4.4.1. Influence of silicate absorption on PAH ratios 

The fluxes of the PAH bands at 6.2/<m, 7.7/im, 8.6/im, 
and 11.3/im are all affected by dust absorption. However, 
as the 11.3/im is significantly more affected by the 9.8/im 
silicate feature than the 6.2/im feature, this will affect the 
measured PAH ratio in Fig. [H To illustrate this influ- 
ence, we show in Fig. [S] the baseline subtracted spectra 
of the selected regions around the nucleus in M82, nor- 
malized to the flux at 7.7/im to emphasize the relative 
flux v ariations of the 11.3/iin PAH feature. The gray line 
is the IChiar fc TielensI (|2006r) extinction law, applied to 
a normalized flux of 1.75, and assuming rg.g = 2. For 
this figure, the baseline was removed by subtracting a 
second order polynomial fitted to the following wave- 
lengths: 5.5/im, 6.8/im, 8.0/im, 13.2/im, and 14.5/im. 
These wavelengths were chosen to avoid the silicate fea- 
ture at 10/im and the PAH features. The effect of ex- 
tinction on the 8.6/im and 11.3/im PAH features is simi- 
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Fig. 9. — Baseline subtracted IRS spectra of the selec ted regions 
from Fig. 1. The gray line is the lChiar Hi TielensI II2006I) extinction 
law, applied to a normalized flux of 1.75, and assuming rg.g = 2. 

lar, meaning that the 8.6/11.3/xm ratio could possibly be 
used to study PAH ionization with minimal concern for 
extinction, as the difference in flux correction between 
8.6/ini and 11.3/im is ~ 15%. However, the 8.6/im fea- 
ture is influenced by the broad PAH feature at 7.7/im, 
which makes the local continuum fitting more difficult 
compared to the 6.2/xm feature. 

Fig. [9] suggests that the variations of PAH feature ra- 
tios involving the 11.3/im feature will be heavily affected 
by extinction and contributing to the distribution of the 
6.2/11.3^m PAH ratio in Fig. [S] The 7.7/zm PAH fea- 
ture is affected by extinction to a similar level as the 
6.2//n i feature, as shown in the Fig. [9l iDraine 
()2001| ) use the 6.2/7.7/im ratio for PAH size diagnostic 
and 11.3/7.7/xm for a PAH ionization state diagnostic. 
Variations of these ratios reflect real variations of the 
physical properties of PAHs in M82 only if their varia- 
tions are not due to extinction effects. Fig. [TUl presents 
the 6.2/7.7/im and 11.3/7.7/im PAH ratios calculated 
from PAHFIT measurements of the same sub-regions as 
in the extinction map. The data in Fig. [TO] is corrected 
for extinction. We took the average of the extinction 
methods (continuum fit method and the PAHFIT fit on 
SL-l-LLl spectra (see 4.3)) for correction. The points are 
dispersed between the tracks representing totally ionized 
and neutral PAH populations, and the error bars show 
the average of the difference between the extinctions de- 
rived by the two methods. The triangles are points from 
the areas where [NeHI] / [NeH] > 0.24 in the [NeIII]/[NeH] 
map from Fig. O These are regions with harder radia- 
tion field, where a greater number of ionized PAH are 
expected to be observed. The squares are points from 
areas in the map where [NeHI]/[NeH]< 0.13. To preserve 
the clarity of the plot, only the error bars at these points 
are represented, as they are typical for all the points. 
The arrow represents the effect on the ratios of a silicate 
absorption feature with rg.g = 1.34, which is the median 
optical depth for the region as explained in section 4.3. 

The 11.3/7.7/xm PAH ratios are dispersed between 0.2 
and 0.5. However, the error bar can account for this dis- 
persion, demonstrating how silicate absorption can affect 
the diagnostic of PAH ionization state based on this ra- 
tio, on scales of the size of the IRS resolution elements. 
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Fig. 10.— The e.2/^.^^J,m PAH ratio vs. 11.3/7.7 PAH ratio, 
for the same sub-regions as in the absorption map in Fig. [7] The 
triangles are points from the areas where [NeIII]/[NoII]> 0.22 and 
the squares are points from the areas where [NeIII]/[NcII]< 0.13. 
Error bars represent average uncertainties, related to extinction 
correction and fitting uncertainties. The arrow represents the ef- 
fect of a silicate absorption feature with rg.g = 1.34. The tracks 
represent the ratios of a population of neutral and ionized PAHs. 

~ 35 pc. This implies that the average variations in the 
ionization state of the grains are relatively small on scales 
of 35 pc. These results are to be compared to SINGS re- 
sults, in which the 11.3/7.7/im PAH ratio is found within 
the same range as in the centr al region of M82 for galax- 
ies dominated by HII regions (jSmith et al.l[2006ah . 

The 6.2/7.7/im ratio varies between ~0.2 and ~0.3. 
This is weff in th e range of 0.2 — 0.4 reported by 
ISmith et al.1 (|2006af ) for SINGS galaxies dominated by 
HII regions. We observe no significant difference in the 
6.2/7.7/im ratio between regions with high [NeIII]/[NeII] 
ratio and regions with less [NeIII]/[NeII] ratio. Extinc- 
tion and fitting errors affect the 6.2/7.7/im ratio only 
by ^ 1%, and errors from fitting residuals from PAH- 
FIT amount to less than 2%. However, there are sig- 
nificant uncertainties arising from the continuum fitting 
by blackbody curves with effective temperatures between 
35 — 300K. Modifying the number of blackbody compo- 
nents and their temperatures results in changes of the 
PAH ratios in excess of 8 — 12%. The horizontal error 
bar represents the average uncertainty of 10%. 

The range of 6.2/7.7/im ratios implies an environment 
composed by a warm ionized medium an d photodissocia- 
tion regions (PDRs) (jPraine fc Lill2001f ). The dispersion 
in the data could reflect a real variation of this ratio. The 
general signi ficance of the variation of the 6.2/7.7/im is 
discussed on Draine & Li (2001). Assumptions about the 
stellar radiation intensity, which affect the 11.3/7.7/im 
ratio as well, account f or these variations in the number 
of C atoms. Following (jPraine lFLil[200l . our observed 
ratios correspond to a number of carbon atoms in a PAH 
grain between 100-140 (6.2/7.7/im- 0.3) and 240-320 
(6.2/7.7/im- 0.2). 

However, there is no correlation between the 6.2/7.7/tm 
ratio with radiation hardness (symbols in Fig. [TO]) or 
with any other resolved spatial structure in the central 
region. Since the variations are comparable to the un- 
certainties in the measurement the results are not (yet) 
signiflcant enough to support PAH size variations at par- 
sec scales. Similarly, the variations of the 11.3/7.7/im ra- 
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tio are mostly due to extinction and show little support 
for variations in PAH ionization throughout the region. 
However, if PAH sizes would vary on parsec scales one 
would expect a strong emission from the larger grains 
which radiate predominantly at larger wavelengths, such 
as the 17 /mum complex. 

4.4.2. The 17^m PAH complex 

The 16 — 18/im wavelength range contains the H2 
8(1)17.03 line and the 17/im PAH complex. This 
complex is attributed to a blend of emission fea- 
tures (16.45^m, 17.03/im, and 17.37/im) which are 
possibly due to PAH C-C -C bending modes (e.g. 
IVan Kerckhoven et al.l l2000( ). and also emission from 
PAH clusters, amorp hous carbon particle s, and other 
PA H-related spe c ies dPeeters et al.1 l2004a| ) . Following 
the iDraine fc Lil ()2007| ) models, the 17^m complex is 
mostly emitted by large PAH molecules with 1000 - 
2000 carbon atoms, while the 6.2/xm feature is emit- 
ted mostly by smaller PAH molecules with only 200 - 
300 carbon atoms. The relatively narrow features at 
16.45/zm and 17.37Mni have been detected by ISO in 
Galactic (iMoutou et al.l[2000l ) and extragalactic sources 
(|Sturm et al.l l2000( l. Only with Spitzer-IRS has the 
whole 17/im PAH complex been identified and routinely 
detected in both normal an d starburst galaxie s (e.g . 
ISmith et al.1 12004 iDale et al.l [2006: Brandl et al. l2006h 
inclu ding the outer regions of M82 (E ngelbracht et al.1 
120061 ). However, the spatial variation within galaxies 
other than the M i lky W ay has not yet been studied. 
lEngelbracht et aT] ()2006l ) detected this complex in the 
disk and in the halo of M82, but their spectral coverage 
did not include the central regions. Here we report the 
characteristics of this complex in the central kpc of M82. 

As seen in Table [3] the ratios between the 17/im com- 
plex and the 11.3/tm PAH feature, as measured with 
PAHFIT, range from 0.39 in regions 1 and 3 to 0.49 in 
the center region, and do not show clear spatial correla- 
tions. However, the 17/6. 2/tm varies from 0.21 in region 
3 to 0.57 in region 4. The ratios between the regions 
in the plane and the outward regions vary by a factor 
of 2. The latter numbers suggest a significant variation 
possibly due to different PAH sizes. However, correcting 
for an optical depth of rg.g = 2 for the regions around 
the [NeH] peaks (Fig. [7|), the ratio increases by ~ 67%, 
yielding values of 0.37 for region 2 and 0.35 for region 3. 
While substantially increased, the values are still more 
than 20% below the values for region 1 and 4. However, 
due to uncertainties on the determination of rg.g we can- 
not diagnose any variation of PAH sizes based on this 
difference. 

Fig. [11] shows a baseline-subtracted spectrum of the 
17/im PAH complex normalized to the average flux be- 
tween 16.2/im and 17.6/fm and rebinned to i? = 300, to 
enhance the broad components. In this figure we can 
see the flux variation in these three main PAH features 
between the regions. The differences in the strength 
of these features are consistent within the uncertainties, 
with exceptions being region 4 at 16.45/im and regions 2 
and 3 at 17.37/im. Moreover, these discrepancies could 
be explained by errors in baseline subtraction, which are 
in the order of 10%. The 17/im PAH complex is also 
seen in the LL spectr a as far as 2 kpc fro m the galactic 
plane. The model of iDraine fc Lil (|2007f) predicts that 
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Fig. 11. — Zoom-in on the PAH dominated emission features in 
the 16 — 18/im spectral range of the SH spectra. The spectra were 
normaUzed to the average flux between 16.2/im and 17.6/tm and 
rebinned to a resolution R = 300, to diminish the noise. The error 
bars indicate noise at three wavelengths. 

the relative strength of the individual components is not 
sensitive to ionization state or grain size. Our observa- 
tions are in good agreement with an invariant shape of 
this PAH complex. 

4.4.3. The difference between PAH emission and VSG 
emission 

Very small grains (VSGs) are dust particles larger on 
average than PAHs, with typic al sizes in the range of 
1 - 150 nm (jPsert et al.lll990l) . They are excited by 
stochastic heating and are thought to be responsible for 
most of the mid-infrared continuum emission. The prop- 
erties of these grains with relation to PAHs have been 
studied previously in galactic sources (jVerstraete et al.l 
1996; Lcbouteillcr et a l] l2007f) and in dwarf galaxies 
f|Madden et al.,.2006 : iWu et al.ll2006[ ). These studies fo- 
cus on the behavior of VSGs in conditions where PAHs 
are destroyed, mainly by intense stellar radiation in HII 
regions or a low metallicity environment. Flux differ- 
ences between VSGs a nd PAHs are reported in Galac- 
tic HII regions by (e.g. iLebouteiller et al] 120071 ). where 
the PAH/VSG emission ratio increases with the distance 
from the cluster. 

As shown in Table 3, there are significant differences 
in EWs of the PAH features between the regions in the 
galactic plane and the regions outside the galactic plane. 
However, the EW shortwards of 10/im behaves differently 
from the EW longward of 10/tm. Longwards of 10/tm, we 
observe that the EWs of the PAH features increase out- 
wards the galactic plane. The 11.3/im EW, for example, 
increases from 0.660/im in region 2 to 2.27/im in region 1. 
The EWs of the 6.2/im and 7.7/im decrease outwards the 
galactic plane. For example, the 6. 2/tm EW decreases 
from 1. 10/tm in region 2 to 0.308/im in region 1. The 
EW of the 8.6/tm feature decreases up to a factor of two, 
from 0.865/im in region 3 to 0.437/tm in region 2. This 
could be due to the contribution of the stellar continuum 
to the local continuum shortwards of 10/tm. The stellar 
continuum contribution decreases with wavelength and 
with the distance from the galactic plane. 

The continuum longwards of 10/tm is composed by 
thermal emission from VSGs. As the PAH flux does 
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not increase outwards from the galactic plane, this can 
only be due to a decrease of the local continuum emis- 
sion relative to the PAH strength. This decrease can be 
explained by several different factors: photo-destruction, 
abundance differences between VSGs and PAHs, and dif- 
ferent heating opacities between VSGs and PAHs with a 
varying radiation field. 

Decreasing PAH flux with the hardening of the ra- 
diation h^s_been_observed_hi galactic star forming re - 
gions (jVerstraete et all Il996t ILebouteiller et all l2007l ). 
on small spatial scales near the luminous clusters (~2 
pc), where the PAH destruction largely surpasses PAH 
excitation. Wether the conditions that lead to PAH de- 
struction on 2 pc scales can be maintained over much 
larger scales, corresponding to the resolution of our maps 
(~35 pc), cannot be derive d from our data, but has 
been observed in NGC 5253 (iBeirao et al.ll2006D. Using 
the SINGS sample of galaxies. iDraine et al. have 
shown that the fraction of PAH abundance over the total 
dust abundance decreases with metallicity. Hence, abun- 
dance differences between PAHs and VSGs could be ob- 
served in cases of a strong metallicity gradient. M82 has 
a metallicity gradient, but it becomes noti ceable only at 
dista nces larger than 1 kpc from the center (jRanalli et al.l 
|2006[) , which is greater than the distance from our regions 
1 and 4 to the galactic plane (400 pc). 

Our favored explanation is the difference in excita- 
tion between PAHs and VSGs, enhanced by variations 
of the radiation field. The mid-infrared PAH features 
are produced by vibrational-rotational modes of the 
PAH molecules, while the VSG continuum emission is 
mainly produced by thermal radiation. VSGs are big- 
ger than PAHs, so thermal radiation becomes domi- 
nant over vi b ration al-rotational transitions. Following 
iDraine fc Lil (j2001l ) models, the opacity of the cross- 
section peaks at FUV wavelengths, where hot dust is 
needed to emit at 10 — 20/im. However, at the distance 
from the galactic plane of regions 1 and 4, the radiation 
field is still intense enough to excite PAHs, but no longer 
of a high enough intensity to excite the dust to the same 
temperatures as in the plane of the galaxy. This provokes 
the PAH EW enhancement observed in these regions. 

4.5. Excitation of the Warm H2 

Molecular hydrogen is the most abundant molecule in 
the Universe and can be used to probe the properties 
of the warm molecular gas in M82. It can be traced in 
the mid- infrared through rotational emission lines, which 
may arise through three different mechanisms: UV exci- 
tation in PDRs surrounding or adjacent to the HII re- 
gions; shocks that accelerate and modify the gas in a 
cloud, coUisionally exciting the H2 molecules; hard X- 
ray photons capable of penetrating the molecular clouds 
and heating large ionizing columns of gas. 

Vibrational-rotational H2 line emission i n M8 2 was 
studied in the near- infrared by iPak et all (j2004D . By 
correlating the emission flux from these lines with [CII] 
157/im and far-infrared luminosity, they showed that the 
H2 emission comes mainly from the PDRs. ISO obser- 
vations of M82 detected S(0), S(l) S(2 ), S(6), and S(7) 
rotational lines (Rigo poulou et ani2002l ). excited by UV 
radiation from massive stars. 

The S(l) and S(2) rotational transition lines of H2 are 
clearly detected in our SH spectra. The fluxes and tem- 
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Fig. 12. — The PAH 11.3^m feature over the 15/im continuum 
versus H2S(1) over the I5fim continuum, based on a pixel-by-pixel 
correlation between SH maps. The values corresponding to each 
of the selected regions are represented by different symbols. The 
error bars represent line flux measurement errors. 

peratures of these lines are listed in Table (H as well as 
the S(l)/S(2) ratios. The temperatures were calculated 
from the S(l)/ S(2) ratios using the method described by 
Rous sel et al.l (j2006). assuming an ortho- to para- ratio of 
three. Our derived temperatures are in agreement with 
the average temper ature value of 450 K, derived with ISO 
(jRigopoulou et al.l f2002). The H2 S(0) line could not be 
detected in any of the LH spectra, as its equivalent width 
is very low, supporting our finding that the H2 temper- 
ature is indeed relatively high. 

If H2 emission is excited mainly by UV radiation in 
PDRs, the fluxes of H2 lines and PAH lines should cor- 
relate closely, since PAH emission features arises from 
the same mechanism in PDRs. To look for secondary 
effects, like shocks, we plot in Fig. |T2]thc PAH 11.3/im 
vs. H2 S(l)17.03/im fluxes, both divided by the contin- 
uum flux at 14.8 — 15.2/im to reduce the range covered 
by the figure. Each data point corresponds to a res- 
olution element of the CUBISM flux maps taken from 
SH spectra. The stars correspond to region 1, squares 
to region 2, diamonds to region 3, and triangles to re- 
gion 4. The error bars represent the line uncertainties 
from the measurement. The correlation in Fig. 1121 shows 
that the excitation mechanisms for both species coincide 
at least on spatial scales of the resolution of the map, 
which is 35 pc. While there is wider dispersion of data 
points corresponding to regions 1 and 4 (the regions asso- 
ciated with outflows) the systematic measurement errors 
account largely for this dispersion. We conclude that 
shocks in the outflow regions have no measurable influ- 
ence on the H2 emission. However, we cannot distinguish 
between UV- and local shock excitation as produced by 
supernovae or energetic outflows, on scales smaller than 
the resolution of our pixels (35 pc). 

5. CONCLUSION 

We presented spatially resolved mid-infrared spectra of 
the central region of M82. The spectra were taken with 
the Spitzer Infrared Spectrograph in both the 5 — 38/im 
low-resolution (i? ~ 65 — 130) and the 10 — 37/im high- 
resolution R ^ 600 module. The high signal-to-noisc and 
the continuous spatial and spectral coverage allowed us 
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to study the nucleus of a starburst galaxy in unsurpassed 
detail. Our goal was to study the physical conditions 
of the interstellar medium and their spatial variations 
within the central kpc of M82. 

Overall, the spectra show the typical features of a star- 
burst: prominent PAH features, silicate absorption, fine- 
structure lines, and a steeply rising continuum. We built 
a spectral map with high-resolution spectra, selected six 
representative regions for spectral extraction, and stud- 
ied the variations of the neon ionic lines and PAH feature 
emission among the regions. 

We attempted to trace the structure of the ionizing 
radiation, and therefore the young stellar population of 
M82 through the diagnostic emission lines, [NeII]12.8/im, 
and [NeIII]15.5/xm. The overall morphology of the ion- 
izing radiation in the region as revealed by the spectral 
maps appears to be in agreement with previou s ground- 
based observations (|Achtermann fc Lacvlll995[ ) . There is 
surprisingly little spatial variation of the [NeIII]/[NeII] 
ratio across the disk, and it varies only by a factor of 
three throughout the central region, with these variations 
not corresponding to the position of the emission peaks. 
We observed an increase of the [NeIII]/[NeII] ratio by a 
factor of 2 outwards the galactic plane, which may be as- 
sociated with the outflows. We suggest that the increase 
of the [NeIII]/[NeII] ratio with distance to the galactic 
plane is due to a decrease in gas density rather than a 
hardening of the field. 

The [Nelll] / [Nell] ratio is low on average, which in- 
dicates that the dominant population consists of already 
evolved clusters (> 5 Myrs). We cannot rule out the pres- 
ence of ongoing star formation in the central region, but 
it must occur at a much reduced rate (< 5%) compared 
to previous epochs. This drop is unlikely to be caused by 
starvation, as there is still a large amount of molecular 
gas present in the central region. It is more likely due 
to negative feedback processes causing a decrease in the 
star formation rate. 

There are significant variations in the amount of dust 
extinction, which strongly correlate with the CO 1-0 
emission across the central region. These variations are 
strong enough to affect the interpretation of PAH fea- 
tures, but due to limited spectral coverage the extinction 
estimates are uncertain. 

The flux of the main PAH features correlates spatially 
with the flux of the neon ionic lines, and with previous 
IRAC observations. Variations in the PAH ratios such 
as 6.2/11.3/xm were observed across the disk. However, 
they are strongly affected by the silicate feature at lO^m. 
We studied the variations the 6.2/7.7/im and 11.3/7.7/im 
PAH ratios, which are diagnostics for the size and the 
degree of ionization of PAHs. We found no correlation 
between the 6.2/7.7/im ratio with radiation hardness or 
with any other resolved spatial structure in the central 
region. Since the variations are comparable to the uncer- 
tainties in the measurement the results are not (yet) sig- 



nificant enough to support PAH size variations at parsec 
scales. Similarly, the variations of the 11.3/7.7/im ratio 
are mostly due to extinction and show little support for 
variations in PAH ionization throughout the region. 

The 17/im PAH complex is very prominent in the cen- 
ter of M82. We did not find any relative variations within 
the complex, which is in agreement with predictions. The 
variations of the 17/6. 2/xm ratio are most likely due to 
extinction effects. Due to the uncertainties on the deter- 
mination of extinction, we did not consider the remaining 
variations as a clear indicator of PAH size variation. 

We observed an enhancement of the EWs of the 
11.3/im and the other PAH features longwards of 10/xm 
outwards from the galactic plane. Several explanations 
exist for this, but we favor the variation of the UV ra- 
diation field, which excites differently PAHs and VSGs, 
given their different sizes. 

The S(l) and S(2) rotational transition lines of H2 have 
been detected in our spectra throughout the central re- 
gion. H2 and PAHs coincide at least on spatial scales of 
the resolution of the map, which is 35 pc. We conclude 
that large scale shocks in the outflow regions have no 
measurable influence on the H2 emission. However, we 
cannot distinguish between UV- and local shock excita- 
tion as produced by supernovae or energetic outflows, on 
scales smaller than the resolution of our pixels (35 pc). 

The Spitzer-IRS observations of the central region of 
M82 complements previous studies, not only in mid- 
infrared, but also in other wavelengths. Our results 
demonstrate the importance of spatially resolved spec- 
troscopy in starburst studies. They helped to constrain 
the age of the starburst and confirm results from other 
studies and also stressed the importance of a thorough 
study of extinction to investigate possible variations of 
PAH properties. Further research of the starburst feed- 
back and quenching processes will elucidate the sharp 
decrease in star formation in the last 5 Myr. It would be 
interesting to see if higher spatial resolution, as expected 
from the MIRI on board of JWST will discover a wider 
variation of the radiation field or signs of variation of 
PAH sizes an/or ionization. In addition, spatially more 
extended spectral studies of M82 are necessary to study 
the connection of the ionic line ratios with the outfiows. 
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TABLE 2 

Main Emission Line Strengths in units of 10~^^W cm~^ 
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1 


0.30±0.03 


0.04±0.01 


1.04±0.03 


0.22±0.01 


0.38±0.01 


0.12±0.01 


0.07±0.01 


1.71±0.46 


381±44 


0.21±0.02 


2 


7.82±0.23 


0.31±0.01 


25.3±0.9 


3.42±0.09 


8.5±0.19 


0.23±0.04 


0.25±0.03 


0.92±0.30 


586±117 


0.13±0.01 


3 


5.77±0.10 


0.33±0.01 


18.4±0.4 


3.58±0.06 


5.92±0.18 


0.25±0.04 


0.27±0.03 


0.93±0.28 


583±109 


0.19±0.02 


4 


0.44±0.02 


0.03±0.01 


1.73±0.04 


0.36±0.01 


0.66±0.02 


0.08±0.003 


0.06±0.003 


1.33±0.13 


439±31 


0.21±0.02 


Center 


17.1±0.47 


0.72±0.08 


55.4±0.97 


8.53±0.15 


18.6±0.37 


0.51±0.07 


0.67±0.07 


0.76±0.21 


700±129 


0.15±0.01 


Total 


44.8±1.49 


2.16±0.13 


141.2±3.10 


23.6±0.38 


49.0±0.83 


2.53±0.09 


2.47±0.11 


1.02±0.09 


536±32 


0.17±0.02 



TABLE 3 

Main PAH Feature Strengths as Measured with PAHFIT 
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0.866 


0.848 


0.660 


0.734 


0.717 


0.753 






3 


136.3±0.1 


478.2±0.2 


83.6±0.1 


73.7±0.1 


57.2±0.1 


5.5±0.1 


29.1±0.1 


0.21 


0.39 




0.827 


0.884 


0.865 


0.784 


0.877 


0.850 


0.720 






4 


13.3±0.1 


49.6±0.1 


11.3±0.1 


15.7±0.1 


7.7±0.1 


0.67±0.10 


7.6±0.1 


0.57 


0.48 




0.390 


0.320 


0.575 


1.86 


1.52 


1.30 


2.16 






Center 


375.0±0.2 


1320.0±0.5 


242.4.0±0.5 


202.5±0.6 


156.3±0.3 


14.3±0.1 


99.0±0.2 


0.26 


0.49 




1.22 


0.991 


0.937 


0.758 


0.791 


0.750 


0.907 






Total 


1158.0±1.0 


4121.7±3.4 


743.9±1.5 


736.0±1.7 


535.8±0.7 


52.2±0.2 


327.8±0.5 


0.28 


0.45 




1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 







* Flux in units of IQ-^^W cm ^. The quoted errors are only the PAHFIT fit residuals. EWs are normalized to the EWs for the whole SH slit.*^ Equivalent 
width in units of /^m 



